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ABSTRACT 


The  objective  of  this  project  is  to  use  computational 
fluid  dynamics  to  examine  the  flow  around  three  Schiebe 
head  forms:  2  inches,  10  inches,  and  20  inches  in  diameter. 
The  flow  around  the  20-inch  head  form  was  examined 
experimentally  in  the  Large  Cavitation  Channel  by  making 
laser  Doppler  velocimetry  measurements.  Computational 
analysis  was  used  to  determine  the  flow  around  all  three 
head  forms.  The  computational  results  and  the  experimental 
results  were  compared.  The  objective  of  the  comparison  was 
to  determine  if  a  computational  model  could  be  used  to 
represent  accurately  a  physical  model  in  a  test  tunnel. 

The  results  showed  that  both  viscous  and  inviscid 
flow  calculations  represented  the  experiment  well.  The 
inviscid  analysis  was  much  less  costly  in  time  and  resources 
while  still  providing  useful  results.  Inviscid  flow  calculations 
seem  to  be  the  best  choice  for  a  preliminary  prediction 
technique  to  support  experimental  investigations. 

This  report  presents  a  description  of  the  experimental 
and  computational  methods  used,  a  detailed  comparison  of 
results,  and  an  analysis  of  the  comparison. 


ADMINISTRATIVE  INFORMATION 

This  work  was  funded  by  the  Large  Cavitation  Channel  project  office  (Code  1502), 
Carderock  Division,  Naval  Surface  Warfare  Center,  headed  by  Dr.  Wm.  Middleton  and  Mr.  P. 
Yamall,  with  Job  Order  Number  1-1502-207-03. 

INTRODUCTION 

The  recent  development  of  the  Large  Cavitation  Channel  (LCC)  has  brought  about  a 
chance  to  implement  state-of-the-art  instrumentation  and  experimental  techniques.  Accurate 
and  intricate  experimental  investigations  of  flow  phenomena  are  made  feasible  by  these 
procedures.  Laser  Doppler  velocimetry  (LDV)  is  among  these  techniques.  Computational  fluid 
dynamics  (CFD)  is  another  state-of-the-art  discipline  that  is  now  leaving  the  "research-only" 
stage  and  entering  practical  use.  Computational  techniques  include  viscous  flow  codes  and 
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inviscid  flow  codes.  Recent  and  on-going  improvements  to  these  algorithms  have  enabled  them 
to  be  used  for  increasingly  accurate  estimates  of  flow  phenomena. 

The  LCC  presents  an  opportunity  to  exploit  the  capabilities  of  these  technologies,  LDV 
and  CFD,  and  use  them  to  complement  each  other.  The  primary  benefits  are  two-fold.  First,  the 
technologies  allow  much  more  detailed  and  complicated  investigations.  Second,  the  traditional 
experimental  method  is  made  more  efficient  and  reliable  by  providing  a  system  for  more  careful 
planning  and  execution  of  measurements. 

The  ultimate  goal  of  this  integration  is  to  develop  a  system  where  computations  are  used 
as  a  preliminary  prediction  technique  for  flow  characteristics.  Physical  investigations  would 
concentrate  in  areas  of  interest  identified  by  the  early  calculations.  A  computational-to- 
experimental  correlation  could  then  be  determined,  and  the  calculations  revised  to  represent  the 
actual  flow  more  accurately.  A  database  of  correlations  would  be  constantly  updated,  leading  to 
increased  accuracy  in  future  preliminary  predictions.  Further  investigations  of  the  same  model, 
including  small  changes  in  model  shape,  appendages,  test  condition,  etc.,  could  be  examined 
computationally,  reducing  the  use  of  costly  experimental  resources  and  facilities. 

Such  a  system  cannot  be  fully  implemented  at  this  time.  The  data  base  to  make  accurate 
and  reliable  correlations  has  not  been  gathered  yet.  Nor  are  the  intricacies  of  the  techniques 
involved  currently  understood  to  the  degree  necessary  to  use  their  potential  fully.  Still,  the 
experience  to  use  these  techniques  must  be  gained  when  possible,  so  that  an  integrated 
computational-experimental  system  can  be  implemented  without  a  long  delay.  The  growing 
costs  in  time  and  money  for  experiments  push  toward  the  use  of  a  more  efficient  system  with 
increasing  power.  This  project  represents  the  first  step  toward  an  integrated  experimental- 
computational  analysis  system. 

MODEL  AND  TEST  FACILITY  DESCRIPTION 

The  three  models  examined  consisted  of  a  Schiebe  head  form  forward  of  an  axisymmetric 
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body  with  e  tapered  tail.  The  Schiebe  series  of  bodies  was  developed  by  calculating  the  potential 
flow  around  a  source  disk  in  a  uniform  field,  yielding  a  minimum  pressure  coefficient  (Cp,min)-  As 
the  source  disk  radius  (B)  is  decreased,  the  minimum  pressure  coefficient  on  the  head  form 
surface  increases,  with  a  point  source  (B=0)  yielding  a  minimum  pressure  coefficient  of  -0.33. 
The  series  of  bodies  is  classified  based  on  minimum  pressure  coefficient  and  geometrically  scaled 
to  the  desired  size.l*  The  three  models  examined  had  geometrically  similar  Schiebe  head  forms 
with  a  minimum  pressure  coefficient  of  -0.75.  The  head  forms  were  2  inches,  10  inches  and 
20  inches  in  diameter.  Table  1  gives  the  non-dimensional  offsets  of  the  Cp,niin=-0-75  Schiebe 
head  form. 

The  head  form  was  followed  by  a  body  of  revolution  in  each  case.  Venical  tail  fins  were 
added  to  the  10-inch  and  20-inch  bodies  to  minimize  bending  moments  on  the  support  strut.  The 
small  cross-section  and  long,  slender  shape  of  the  2-inch  model  was  stable  without  fins.  All 
three  models  were  mounted  on  the  same  ogive  strut.  Figure  1  shows  profile  drawings  of  these 
models. 

The  experimental  analysis  was  undertaken  in  the  test  tunnel  at  the  Large  Cavitation 
Channel  (LCC),  a  Carderock  Division,  Naval  Surface  Warfare  Center  (NSWC)  detachment  in 
Memphis,  Tennessee.  The  test  section  in  the  tunnel  is  40  feet  long,  and  10  feet  by  10  feet  in 
cross-section,  with  1.25  foot  notches  in  each  comer.  Windows  are  located  on  each  surface 
around  the  tunnel  test  section  throughout  the  length.  The  windows  on  the  top  and  bottom  are 
obstructed  by  ceiling  plates  and  an  acoustic  trough,  respectively.  Excellent  visibility  is  afforded 
to  the  tunnel  centerline  along  the  entire  length  from  the  test  platform.  The  structure  around  the 
windows  limits  the  accessibility  of  an  orthogonal  line-of-sight  to  the  model,  which  is  necessary 
when  making  laser  Doppler  velocimetry  measurements. 

In  the  computational  analysis,  the  bodies  were  arranged  such  that  the  model  axis  was 


*  References  are  listed  on  page  61. 
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incident  with  the  tunnel  centerline.  In  the  physical  experiment  of  the  20-inch  model,  the 
longitudinal  axis  was  aligned  about  one-half  degree  nose  down,  and  one-half  degree  yawed 
away  from  the  control  room,  but  very  close  to  the  tunnel  centerline. 

COMPUTATIONAL  ANALYSIS 

The  computational  analysis  process  involves  many  steps.  The  numerical  models  must  be 
developed,  much  like  physicals  models  are  developed.  A  numerical  representation  must  be 
generated,  suitable  to  the  particular  type  of  calculations  being  performed.  The  calculations  are 
made,  and  the  solution  must  be  iterated  until  reasonable  results  are  achieved.  Finally,  the 
results  must  be  presented  in  a  meaningful  form.  This  process  involves  many  computer  programs 
and  tools.  Figure  2  shows  an  overview  of  this  process.  The  basics  of  this  process  are  described 
below.  A  detailed  description  of  the  procedure  appears  in  Appendix  A. 

GEOMETRIC  MODELLING 

To  begin  the  computational  analysis,  numerical  models  were  developed.  Representations 
of  the  cavitation  head  form  models,  the  strut,  and  the  LCC  test  section  were  generated  from 
assembly  drawings  and  offsets.  Three-dimensional  numerical  surfaces  representing  the  desired 
physical  model  were  generated.  As  the  surfaces  were  defined,  a  grid  was  distributed  over  the 
surface.  Shaded  pictures  of  these  surfaces  were  generated,  as  shown  in  Figure  3.  These 
pictures  are  commonly  used  for  presentations,  discussions,  and  visual  inspections. 

As  is  common  for  symmetrical  forms,  only  half-models  were  developed  for  the 
computations.  For  the  viscous  flow  analysis,  a  model  of  the  20-inch  test  body  was  generated 
with  the  strut  and  tunnel  walls.  For  the  inviscid  flow  computations,  a  model  of  the  20-inch  body 
and  strut  was  generated  both  with  and  without  the  tunnel  walls.  The  tunnel  walls  did  not 
influence  the  flow  noticeably  in  the  primary  area  of  interest  (i.e.,  the  head  form),  as  explained 
later.  Therefore,  the  inviscid  panel  models  for  the  2-inch  and  10-inch  bodies  were  generated  with 
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the  strut  but  without  the  tunnel  walls. 

VISCOUS  FLOW  CALCULATIONS 

Grid  Ggneration 

A  three-dimensional,  multi-block,  volume  grid  was  developed  around  the  surface  grid.  A 
block  structure  was  developed  around  the  desired  geometry.  A  sample  of  the  block  structure  is 
shown  in  Figure  4.  The  heavy  lines  in  the  figure  represent  block  edges.  The  half-model  is  clearly 
visible,  as  is  the  outline  of  the  tunnel. 

The  topology  of  the  blocks  was  defined,  as  was  the  connectivity  of  the  blocks.  Figure  5 
shows  selected  block  faces  from  the  grid  developed  for  the  20-inch  model  in  the  tunnel.  The 
heavy  lines  in  the  figure  again  represent  block  edges,  while  the  thin  lines  show  the  grids  on  the 
block  faces.  For  clarity,  only  a  few  block  faces  are  shown. 

Once  the  face  grids  were  defined,  the  coordinates  of  the  points  on  the  interior  of  each 
block  were  interpolated.  These  points  outline  cells,  or  three-dimensional  volumes.  The  cells 
were  checked  to  ensure  that  their  orientation  was  consistent  (i.e.,  in  a  right-handed  coordinate 
system).  Iterations  of  the  grid  generation  were  made  to  achieve  an  acceptable  three- 
dimensional  volume  grid  with  a  relatively  smooth  and  properly  oriented  grid  distribution. 

The  details  for  developing  a  block  structure  and  grid  distribution  vary  greatly  from  case  to 
case,  and  are  sometimes  based  heavily  on  trial -and-error.  Each  model  has  different  properties 
which  define  various  characteristics  of  the  block  structure  and  grid  distribution.  Hopefully,  early 
determination  these  characteristics  can  be  made  easier  in  the  future.  For  practical  application, 
criteria  must  be  developed  to  help  guide  the  grid  generation.  The  riteria  will  be  based  on  the 
collective  and  continuing  experience  from  successfully  developed  grids.  In  addition,  the  constant 
refinement  of  existing  software,  along  with  the  development  of  new  software  and  hardware, 
facilitates  the  process  of  numerical  model  generation. 

The  20-inch  body  was  modelled  with  the  vertical  tail  fins  and  mounted  on  the  strut.  The 
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inner  boundary  of  the  volume  grid  was  the  body  surface  and  tunnel  centerline,  and  the  outer 
boundaries  were  the  tunnel  walls.  For  computational  reasons,  the  grid  was  extended  forward  of 
the  model  by  one-half  a  body  length,  and  aft  by  a  full  body  length.  Figure  4  shows  the  grid 
outline  well. 

Solution  Technique 

A  steady-state,  incompressible  viscous  flow  calculation  was  made,  using  the  three- 
dimensional  Reynolds-averaged  Navier-Stokes  (RANS)  equations. 

Since  single  grid  can  require  a  total  of  250  to  350  hours  of  computer  time  to  calculate  a 
final  solution,  viscous  flow  analyses  are  potentially  expensive  in  computer  time.  Most  of  this 
time  is  spent  in  the  iterative  process  of  refining  the  grid  for  an  improved  flow  solution.  The  cost 
of  computer  time  and  equipment  corresponds  to  the  material  and  facilities  costs  in  physical  tests. 
Viscous  flow  calculations  currently  can  be  competitive  in  cost  to  physical  model  tests.  The  costs 
of  computational  analysis  become  more  favorable  as  faster  and  less  expensive  computers 
become  available. 

Examining  changes  in  model  configuration  is  often  less  costly  when  done  computationally. 
As  long  as  the  block  topology  and  grid  distribution  does  not  need  to  change  radically,  subsequent 
model  configurations  and  conditions  often  take  only  100  to  150  hours  total  computer  time. 
Physical  model  experiments  generally  require  as  much  time  for  each  configuration,  which  can 
monopolize  facilities,  further  increasing  costs.  parametric  study  of  a  hull  form  or  appendage 
suite  can  be  very  costly  when  tested  physically.  The  same  study  can  be  much  more  extensive 
and  affordable  when  performed  computationally.  Preliminary  and/or  follow-up  physical  tests  can 
be  used  to  verify  computational  results,  thereby  ensuring  the  accuracy  of  the  analysis. 
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INVISCID  FLOW  CALCULATIONS 

Model  Generation 

Inviscid  flow  calculations  are  typically  made  for  one  or  more  panelled  surfaces.  The 
panelled  surfaces  are  generally  the  surface  grids  developed  earlier,  but  with  a  refined  distribution. 
Figure  6  shows  the  panel  distribution  used  for  the  20-inch  model  and  strut.  Note  the 
concentration  of  panels  about  the  area  of  high  curvature  in  the  head  form. 

The  20-inch  model  was  generated  both  with  and  without  the  tunnel  walls.  After  viewing 
the  results,  it  was  determined  that  the  tunnel  walls  did  not  affect  the  flow  over  the  head  form  in 
the  inviscid  flow  solution.  The  tunnel  walls  only  influenced  the  flow  along  the  strut  and  near  the 
model-strut  interface.  The  tunnel  walls  were  not  included  in  the  panel  models  for  the  2-inch  and 
10-inch  head  forms  to  reduce  computation  time. 

Solution  Technique 

The  inviscid  flow  calculations  were  made  for  the  head  form  models.  An  inviscid  flow 
solution  can  take  anywhere  from  2  to  5  minutes  of  computational  time.  Even  multiple  iterations 
while  refining  the  panel  distribution  do  not  involve  much  computer  time.  Viscous  flow  effects, 
such  as  turbulence  and  separation,  are  not  accounted  for  by  inviscid  flow  analyses.  These  effects 
are  often  very  significant,  and  cannot  always  go  unexamined.  On  the  other  hand,  inviscid  flow 
programs  have  the  advantage  of  including  free-surface  effects  more  easily. 

Some  CFD  codes  are  hybrids:  the  program  uses  different  solution  techniques  on  different 
sections  of  the  model.  This  arrangement  has  both  advantages  and  disadvantages.  No  programs 
of  this  type  were  evaluated  during  the  course  of  this  investigation.  They  are  mentioned  only  for 
completeness. 

PRESENTATION  OF  RESULTS 

The  way  in  which  computational  results  are  presented  can  be  nearly  as  important  as  the 
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calculations  themselves.  Tools  are  available  to  display  the  results  in  many  formats.  One  formal 
is  three-dimensional  color  or  shaded  pictures.  Two  figures  show  results  from  the  viscous  flow 
calculations.  Figure  7  shows  a  plot  of  pressure  coefficients  on  the  20-inch  model  and  strut,  while 
Figure  8  shows  particle  traces  around  the  same  model.  For  the  particle  trace,  three  particles 
were  released  near  the  head  form  surface.  Their  traces  were  computed  forward  and  backward  in 
time,  showing  where  the  particles  were  headed,  as  well  as  where  they  had  been.  Three  particles 
were  also  released  just  forward  of  the  upper  fin  leading  edge.  Their  traces  were  computed  aft 
only. 

Results  were  also  plotted  in  two-dimensional  formats  using  various  commercially 
available  graphing/plotting  packages.  Three-dimensional  pictures  are  often  not  the  most  useful 
method  of  presenting  results  for  computations.  Two-dimensional  graphs  are  frequently  used  in 
cases  where  specific  quantities  are  more  important  than  general  trends,  or  cases  where  specific 
comparisons  must  be  made,  as  was  the  case  in  comparing  the  velocity  profiles  and  pressure 
coefficients  in  this  report. 


EXPERIMENTAL  ANALYSIS 

SETUP 

The  laser  Doppler  velocimetry  (LDV)  measurements  were  conducted  only  for  the  20-inch 
model  due  to  time  conitraints  on  the  Large  Cavitation  Channel  (LCC)  facility. 

The  locations  of  test  points  for  the  LDV  system  are  limited  by  tunnel  structure,  lens  focal 
length,  and  the  limits  of  the  traverse  mechanism.  In  this  case,  the  test  points  were  restricted  to 
the  volume  around  the  head  form.  The  physical  extents  of  the  traverse  mechanism  restricted 
movement  past  the  pair  of  windows  at  the  forward  end  of  the  model.  No  points  could  be  located 
immediately  forward  of  the  stagnation  point  due  to  obsU'ucting  window  structure.  The  focal 
length  of  the  lens  allowed  test  points  to  be  located  slightly  past  the  centerplane  of  the  tunnel. 
Figure  9  describes  these  limits  graphically. 
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A  test  matrix  was  developed  to  examine  several  areas  around  the  head  form  at  three 
different  speeds.  The  matrix  allowed  for  testing  above  and  below  the  head  form  in  the 
centerplane,  as  well  as  around  the  model  girth.  A  forward  profile  and  a  profile  from  the  model  to 
the  near  tunnel  wall  were  also  planned.  About  400  samples  were  taken  at  each  test  point. 
Table  2  describes  the  experimental  test  matrix. 

TEST  PROCEDURE 

For  this  particular  LDV  setup,  it  was  determined  that  approximately  350  samples  at  each 
point  must  be  collected  before  the  mean  velocity  becomes  a  steady  value,  as  shown  in  Figure  10. 
The  two  groups  of  samples  shown  in  the  figure  represent  one  highly  turbulent  point  (near  the 
tunnel  wall)  and  one  point  with  low  turbulence  (near  the  model  surface).*  Note  that  even  after 
only  100  samples,  the  mean  is  quite  steady  (within  one  standard  deviation  of  the  mean  velocity). 

The  system  installed  at  the  LCC  was  used  in  a  one-component  back-scatter  mode  for 
this  test.  One-component  mode  indicates  that  only  the  longitudinal  component  of  velocity  (Vx) 
was  measured.  The  receiving  optics  used  were  located  in  the  same  housing  as  the  transmitting 
optics,  so  that  only  light  bouncing  back  toward  the  lens,  or  back-scatter,  was  received.  The 
transmitting  lens  has  a  fixed  focal  length,  and  was  mounted  on  a  three-axis  traverse.  The 
traverse  allowed  the  lens  to  be  moved  in  any  direction,  and  therefore  the  test  point  location  could 
be  moved  in  three  dimensions. 

A  discussion  of  LDV  theory  and  operation  is  beyond  the  scope  of  this  report,  but  is 
readily  available  from  many  sources,  including  Durst  et  alP- 


*  Turbulence  in  any  direction  is  defined  as  the  standard  deviation  of  the  velocity  in  that 
direction  divided  by  the  mean  longitudinal  velocity,  or 

%turbulence  =  =  xl00%  . 

Vx 
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DATA  ANALYSIS 


The  data  were  analyzed  by  averaging  the  valid  samples  for  each  point  (400  samples  in 
most  cases).  Samples  outside  four  standard  deviations  from  the  mean  were  not  included.  This 
filtering  removed  outliers*  without  affecting  samples  in  the  expected  distribution.  Figure  1 1 
shows  examples  of  the  distribution  of  samples  as  well  as  showing  the  filter  criteria.  A  high  and  a 
low  turbulence  point  are  shown.  The  distributions  of  other  points  were  similar. 

The  resulting  mean  velocities  were  graphed  as  velocity  profiles.  Three  velocity  profiles 
are  shown  in  Figure  12.  Tabulated  results  for  all  profiles  are  shown  in  Table  B-1  of  Appendix  B. 


RESULTS  AND  DISCUSSION 

SCHIEBE  RESULTS 

A  potential  flow  analysis  was  used  by  Schiebe  to  determine  the  pressure  distribution  on 
the  body  surface.  ^  Since  the  body  is  axisymmetric,  the  pressure  coefficients  vary  only  with 
longitudinal  distance  from  the  stagnation  point.  The  results  from  Schiebe  are  shown  in  Table  3. 


INVISCID  FLOW 

Inviscid  Flow  versus  Schiebe 

Figure  13  shows  pressure  coefficient  versus  longitudinal  distance  from  the  stagnation 
point  for  the  20-inch  head  form  model,  along  with  the  Schiebe  results.  These  pressure 
coefficients  are  also  listed  in  Table  4.  The  coefficients  were  computed  along  the  side  of  the 
models  at  90  degrees  from  the  upper  edge  and  at  the  lower  edge  along  the  centerline.  Since  the 
side  of  the  model  is  well  removed  from  either  the  strut  or  the  centerline  (a  mathematical 
boundary),  the  results  at  the  side  would  be  expected  to  best  match  the  results  from  Schiebe 
(which  are  for  an  unappended  head  form).  The  comparison  is  very  favorable  for  both  the  side  and 


*  Outliers  are  "readings  that  appear  to  be  significantly  out  of  line  with  the  other  readings."^ 
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the  lower  edge,  however.  Figure  13a  shows  the  pressure  coefficients  along  the  surface  of  the 
head  form,  while  Figure  13b  extends  the  graph  to  show  the  length  of  the  model.  The  results  from 
Schiebe  were  calculated  for  the  head  form  only  (i.e.,  no  strut  or  fins.)  The  results  along  the  lower 
edge  are  not  relevant  to  the  comparison  in  way  of  the  tail  fins,  since  there  is  no  fin  at  the  side  of 
the  model. 

Figure  14  compares  the  pressure  coefficients  of  the  20-inch  half-model  for  the  two 
configurations  tested:  with  and  without  the  tunnel  walls.  Figure  14a  shows  that  there  is  no 
apparent  wall  effect  at  the  side  of  the  model,  while  Figure  14b  shows  the  pressure  coefficients  on 
the  strut  centerline.  In  the  first  graph,  the  wall  effect  can  be  seen  as  a  slight  offset  in  way  of  the 
strut.  The  second  graph  shows  a  difference  in  pressure  coefficient  of  about  0.04  at  the  bottom  of 
the  strut  centerline.  These  data  are  also  presented  in  Table  5. 

This  slight  effect  on  the  model  is  important  to  consider  because  the  wall  definition 
requires  about  the  same  number  of  panels  as  the  model  and  strut  combined,  and  therefore  the 
solution  takes  about  twice  as  much  computer  time.  The  small  effect  at  the  strut  does  not  seem  to 
influence  the  results  on  the  head  form  appreciably.  Since  the  20-inch  head  form  will  have  the 
greatest  wall  effect  of  the  three  head  form  models,  it  is  reasonable  to  assume  that  no  walls  need 
to  be  modelled  for  the  2-inch  and  10-inch  head  forms. 

This  minor  wall  effect  must  be  taken  in  context,  however.  This  axisymmetric  model  was 
20  inches  in  diameter  in  a  tunnel  10  feet  by  10  feet  in  cross-section.  The  ratio  of  cross-sectional 
areas  (i.e.,  blockage  factor)  is  only  '744.  In  addition,  the  model  was  mounted  so  that  its  surface 
was  slightly  over  five  times  the  head  form  radius  from  the  tunnel  walls.  A  test  of  wall  effect 
should  be  carried  out  for  models  of  equal  or  greater  blockage  factor  or  for  models  mounted  closer 
to  walls.  Of  course,  surface  ship  models  mounted  at  the  ceiling  plate  should  always  be  modelled 
with  tunnel  walls. 
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The  wall  effect  is  significant  on  the  strut.  Therefore,  calculations  for  models  in  which  the 
flow  on  or  around  the  strut  are  important  should  include  walls.  Likewise,  the  strut  effect  is 
significant  on  the  model.  Mounting  hardware  and  objects  near  or  on  the  body  can  significantly 
influence  the  flow  characteristics  and  should  always  be  modelled. 

Boundary  Condition  Effects 

Figure  15  presents  the  half-model  results  at  the  lower  edge  of  the  model  and 
corresponding  results  from  a  bare  (unappended)  20-inch  full-model.  The  numerical  definition  of 
the  full-model  starts  and  ends  along  the  upper  edge  of  the  model,  meaning  that  the  lower  edge  is 
180  degrees  removed  from  any  mathematical  boundary.  The  mathematical  boundary  condition 
seems  to  have  little  effect  on  the  solution. 

Figure  13  can  also  be  used  to  examine  boundary  condition  effects.  The  results  at  the  side 
of  the  model  (90  degrees  removed  from  a  boundary)  and  the  lower  edge  (at  a  boundary)  agree 
closely.  The  upper  edge  of  the  model  was  not  considered  due  to  the  presence  of  the  strut. 

This  consistency  is  an  indication  that  the  mathematical  boundary  condition  does  not 
adversely  affect  the  inviscid  flow  calculation  results.  The  full-model  has  twice  the  number  of 
panels  and  requires  about  twice  the  computer  time.  The  half-model  is  a  satisfactory 
representation.  Unless  otherwise  stated,  references  elsewhere  in  this  report  to  the  model  refer 
to  the  half-model. 

Scalc-Effgcts 

Figure  16  shows  the  pressure  coefficients  for  each  of  the  head  forms  at  the  side  of  each 
model.  The  results  for  all  three  head  forms  agree  very  well  with  each  other,  as  well  as  with  the 
Schiebe  results.  Inviscid  flow  theory  deals  primarily  in  non-dimensional  values,  and  would  not 
be  expected  to  show  any  discrepancy  between  the  head  forms.  These  results  are  tabulated  in 
Table  6. 
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It  should  be  observed  that  the  panels  for  the  three  head  forms  were  not  distributed 
similarly,  as  shown  by  the  differing  locations  of  the  data  points  in  the  figure.  The  inviscid  flow 
calculation  results  do  not  vary  significantly  with  small  differences  in  panel  distribution.  The 
ability  to  achieve  consistent  results  when  varying  panel  distributions  is  important  for  a 
preliminary  prediction  technique. 

VISCOUS  FLOW 

Viscous  Flow  versus  Schiebe 

Figure  17  shows  the  pressure  coefficients  for  the  viscous  flow  calculation,  plotted  with  the 
inviscid  flow  calculation  results  and  results  from  Schiebe.  The  viscous  flow  calculation  results 
are  also  listed  in  Table  7.  The  pressure  coefficients  are  compared  at  the  side  of  the  20-inch  head 
form.  The  results  of  the  viscous  flow  calculation  are  not  as  close  to  the  Schiebe  results  as  are 
the  inviscid  flow  calculation  results.  However,  the  results  show  the  correct  trend  and  the 
magnitudes  seem  reasonable. 

A  basic  difficulty  in  working  with  viscous  flow  calculations  is  that  the  results  are  often 
dependant  on  grid  distribution.  This  dependence  inhibits  the  usefulness  of  viscous  flow 
calculations  as  a  preliminary  prediction  technique.  Without  a  prior  estimate,  it  is  difficult  to 
determine  whether  the  viscous  flow  calculation  results  are  reasonable  indicators  of  the  actual 
flow,  or  merely  the  results  of  a  particular  grid  distribution.  In  this  case,  it  would  seem  that  further 
refinement  of  the  grid  in  the  forward  area  of  the  head  form  would  be  important  if  the  grid  results 
were  to  be  used  for  other  investigations.  One  example  of  such  investigations  would  be  a 
parametric  study  to  optimize  the  pressure  disuibution  on  the  head  form.  The  results  are  close 
enough  to  indicate  that  reiteration  of  the  grid  might  yield  more  accurate  results,  and  the  scope  of 
this  investigation  does  not  warrant  further  refinement. 
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Boundary  Condition  Effects 


Figure  18  shows  pressure  coefficient  versus  longitudinal  distance  from  the  stagnation 
point  for  viscous  flow  calculations,  similar  to  the  investigation  presented  in  Figure  13  for  inviscid 
flow  calculations.  The  results  at  the  side  and  at  the  lower  edge  are  compared.  The  comparison 
for  the  viscous  flow  calculation  results  are  not  as  favorable.  The  lower  edge  boundary  condition 
seems  not  to  represent  the  symmetric  plane  very  well  in  the  head  form  area.  Past  the  head  form 
section  (i.e.,  x/L  >  5),  the  side  and  lower  edge  results  agree  quite  well.  The  conclusion  is  that  if 
the  area  near  a  mathematical  boundary  is  of  particular  interest,  full-body  models  should  be 
considered.  The  trade-off  is  computational  cost  and  speed  versus  accuracy. 

The  full-body  model  was  not  tested  for  the  viscous  flow  calculations  due  to  time  and  cost 
considerations.  Likewise,  wall  effects  and  scale  effects  were  not  examined  for  viscous  flow 
calculations. 

Viscous  Flow  versus  Experiment 

Figure  19  shows  a  comparison  of  experimental  results  and  viscous  flow  calculation 
results.  The  results  are  similar  in  trend  and  magnitude  except  near  the  model  surface  and  tunnel 
walls.  The  agreement  is  quite  reasonable. 

In  each  of  the  profiles  compared,  it  appears  that  the  boundary  layer  growth  in  the  viscous 
flow  calculation  is  too  large.  The  turbulence  modelling  in  the  viscous  flow  calculations  uses  an 
empirical  correlation.  The  representation  can  be  changed,  although  sufficient  experimental  data 
must  be  available  to  guide  the  change.  The  results  of  these  investigations  and  future 
investigations  can  be  used  to  compile  a  data  base  for  increasing  the  accuracy  of  the  empirical 
correlation. 

The  modification  of  empirical  correlations  is  not  trivial,  but  could  be  used  if  closer 
agreement  between  the  numerical  and  the  physical  model  was  necessary.  The  cost  of  developing 
a  data  base  and  making  a  modification  was  not  justified  for  this  investigation.  The  cost  would  be 
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justified  for  long  term  computational  support  of  physical  experiments.  Also,  as  more  of  a  data 
base  is  gathered  over  time,  the  cost  of  applying  a  modification  of  this  type  decreased,  while  the 
improved  accuracy  increases,  both  due  to  a  larger  base  or  experience  to  guide  the  change. 

CONCLUSIONS  AND  RECOMMENDATIONS 

The  results  of  this  project  indicate  that  accurate  representations  of  physical  models  can  be 
achieved  with  currently  available  computational  tools  and  techniques.  Computational  fluid 
dynamics  can  reasonably  support  physical  experiments,  and  will  become  more  useful  in  that  role 
as  further  experience  is  gained. 

The  inviscid  flow  calculations  provide  meaningful,  reasonable  data  at  a  relatively  low  cost 
in  time  and  resources.  The  calculations  are  stable  with  regard  to  panel  distribution.  These 
characteristics  enable  the  inviscid  flow  model  to  be  a  competent  preliminary  prediction  technique 
for  physical  experiments  in  a  tunnel.*  Inviscid  flow  calculations  can  provide  the  engineer  with 
vital  information  prior  to  the  experiment. 

For  many  cases,  a  half-model  representation  with  mounting  hardware  is  sufficient  for 
inviscid  flow  calculations.  The  tunnel  walls  may  not  have  to  be  modelled.  Before  these 
assumptions  are  made,  however,  areas  of  interest  should  be  identified  and  the  effects  of  these 
assumptions  considered.  These  assumptions  have  the  potential  to  reduce  computational  time, 
but  could  adversely  affect  the  results  in  certain  situations. 

The  viscous  flow  calculations  can  provide  different  information  than  inviscid  flow 
calculations,  but  require  more  time  and  resources.  Results  are  sensitive  to  grid  distribution, 
making  viscous  flow  calculations  poorly  suited  to  be  used  as  a  preliminary  prediction  technique  at 
the  current  time.  The  development  of  grid  distribution  criteria  is  crucial  to  the  practical 


*  No  effort  was  made  to  validate  free  surface  modelling  capabilities  of  the  inviscid  flow 
calculations,  and  therefore  no  conclusions  can  be  made  as  to  the  proficiency  of  predicting 
experimental  flow  in  non-tunnel  environments  (e.g.,  towing  tanks). 
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application  of  viscous  flow  calculations.  Experience  and  improvements  in  the  techniques  will 
certainly  make  a  difference  in  the  future.  The  results  from  this  investigation  add  to  the 
experience  base  and  can  be  used  for  guidance  in  future  work. 

Currently,  viscous  flow  calculations  seem  best  suited  to  parametric  studies.  A  baseline 
grid  distribution  can  developed  and  results  verified  with  experiments  or  other  prediction 
techniques.  Parametric  changes  could  then  be  made  to  the  grid  without  requiring  the  block 
structure  to  be  significantly  altered.  Also,  cases  where  inviscid  flow  calculations  wilt  not  yield 
the  necessary  information  (e.g.,  separated  flow  characteristics)  are  candidates  for  examination 
using  viscous  flow  calculations. 

Experimental  results  are  important  to  help  guide  the  empirical  techniques  used  in 
turbulence  modelling  for  viscous  flow  computations.  For  example,  experimental  data  should  be 
gathered  to  determine  boundary  layer  characteristics  and  growth.  Once  enough  information  is 
available,  the  empirical  turbulence  model  currently  used  can  be  updated  to  more  accurately 
represent  the  characteristics  of  the  boundary  layer  in  a  tunnel.  Correlation  with  experiments 
would  be  improved,  increasing  the  usefulness  of  viscous  flow  calculations  for  experimental 
support. 

Based  on  these  conclusions,  the  authors  make  the  following  recommendations: 

•  Inviscid  flow  calculations  should  be  used  to  develop  preliminary  predictions  of  flow 
characteristics  prior  to  performing  experiments.  The  cost  of  the  calculations  is 
offset  by  the  availability  of  preliminary  estimates  prior  to  the  experiment.  The 
reliability  of  these  calculations  seems  to  justify  their  use. 

•  Whenever  time  and  resources  permit,  viscous  flow  calculations  should  be 
performed  in  conjunction  with  the  experiment  to  aid  in  the  development  of  tools  and 
techniques.  The  experience  will  help  develop  grid  generation  guidelines,  and  the 
correlation  with  the  experimental  results  will  help  refine  empirical  techniques  used 
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in  viscous  flow  calculations.  These  developments  will  make  viscous  flow 
calculations  less  costly  and  more  reliable  in  the  future  for  experimental  and  design 
support. 

•  The  integration  of  computational  fluid  dynamics  and  physical  experimental 
techniques  should  be  continued  in  order  to  improve  the  overall  experimental 
technique.  This  is  crucial  to  the  modernization  of  experimentation. 
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Table  1.  Schiebe  head  form  offsets  for  the  Cp,min=-0  ’75  head  form. 


Longitudinal  distance 
from  stagnation 
point  (x/R) 

Radius 

(r/R) 

0.00000 

0.00719 

0.00000 

0.01437 

0.00000 

0.02156 

0.00000 

0.02875 

0.00000 

0.03593 

0.00252 

0.14853 

0.00503 

0.20863 

0.00754 

0.25373 

0.01005 

0.29098 

0.01257 

0.32307 

0.02513 

0.44066 

0.03769 

0.51915 

0.05025 

0.57508 

0.06281 

0.61532 

0.07537 

0.64424 

0.08793 

0.66549 

0.10049 

0.68196 

0.11305 

0.69549 

0.12561 

0.70711 

0.14561 

0.72299 

0.16561 

0.73666 

0.18561 

0.74876 

0.20561 

0.75968 

0.22561 

0.76964 

0.24561 

0.77882 

0.26561 

0.78733 

Longitudinal  distance 
from  stagnation 
point  (x/R) 

Radius 

(r/R) 

0.28561 

0.79527 

0.30561 

0.80271 

0.32561 

0.80970 

0.34561 

0.81630 

0.36561 

0.82253 

0.38561 

0.82844 

0.40561 

0.83405 

0.42561 

0.83939 

0.52561 

0.86263 

0.62561 

0.88134 

0.72561 

0.89671 

0.82561 

0.90948 

0.92561 

0.92021 

1.02561 

0.92930 

1.12561 

0.93704 

1.52561 

0.95870 

1.92561 

0.97131 

2.32561 

0.97913 

2.72561 

0.98424 

3.12561 

0.98773 

3.52561 

0.99020 

3.92561 

0.99201 

4.32561 

0.99337 

4.72561 

0.99441 

5.12561 

0.99523 

All  values  are  non-dimensional  by  maximum  head  form  radius. 
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Table  2.  Experimental  test  matrix. 


Profile 

Along  lower  edge. 

forward  of  Cp  min 

4  in.  aft 

8  in.  aft 

16  in.  aft 

Forward  of  stagnation  point 

Around  girth, 

8  in.  aft 

at  upper  edge 
at  45°  above  CL 
at  CL 

at  45°  below  CL 
at  lower  edge 

Along  upper  edge. 

forward  of  Cp,min 

4  in.  aft 

8  in.  aft 

16  in.  aft 

At  CL,  from  model  to  tunnel  wall 

Nominal  tunnel  velocity 


13-6  knots  22.4  knots  29.2  knots 


All  distances  referenced  to  stagnation  point. 
Angles  are  referenced  from  horizontal. 
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Table  3.  Pressure  coefficienis  on  Schiebe  hf^aci  form  (Cp,inin=-0-75)  from  potential  flow  theory. 


Longitudinal  distance 
from  stagnation 
point  (x/R) 

Pressure  coctficient 
(Cp) 

0.00000 

0.99890 

0.00000 

0.99930 

0.00000 

0.99960 

0.00000 

0.99582 

0.00000 

0.99996 

0.00252 

0.98048 

0.00503 

0.95995 

0.00754 

0.93837 

0.01005 

0.91560 

0.01257 

0.89158 

0.02513 

0.74947 

0.03769 

0.56043 

0.05025 

0.31004 

0.06281 

-0.00163 

0.07537 

-0.32280 

0.08793 

-0.56007 

0.10049 

-0.68488 

0.11305 

-0.73585 

0.12561 

-0.75001 

0.14561 

-0.74129 

0.16561 

-0.71918 

0.18561 

-0.69385 

0.20561 

-0.66857 

0.22561 

-0.64444 

0.24561 

-0.62174 

0.26561 

-0.60050 

Longitudinal  distance 
from  stagnation 
point  (x/R) 

Pressure  coefficient 

(Cp) 

0.28561 

-0.58062 

0.30561 

-0.56198 

0.32561 

-0.54447 

0.34561 

-0.52797 

0.36561 

-0.51238 

0.38561 

-0.49761 

0.40561 

-0.48359 

0.42561 

-0.47024 

0.52561 

-0.41187 

0.62561 

-0.36411 

0.72561 

-0.32402 

0.82561 

-0.28983 

0.92561 

-0.26035 

1.02561 

-0.23472 

1.12561 

-0.21232 

1.52561 

-0.14631 

1.92561 

-0.10503 

2.32561 

-0.07811 

2.72561 

-0.05989 

3.12561 

-0.04713 

3.52561 

-0.03793 

3.92.561 

-0.03111 

4.32561 

-0.02593 

4.72561 

-0.02192 

5.12561 

-0.01876 

Distances  are  non-dimensional  by  maximum  head  form  radius. 
Pressure  coefficients  are  referenced  to  the  uniform  flow. 
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Table  4.  Inviscid  flow  calculation  results  for  the  20-inch  head  form. 


Longitudinal 
distance  from 
stagnation  point 
(x/R) 

Pressure  coefficient  (Cp)  | 

At  side 

At  lower 
edge 

0.000 

WBBM 

0.000 

0.000 

0.000 

mmSM 

0.001 

0.9930 

0.001 

mBM 

0.003 

■BH 

0.004 

0.9694 

0.006 

0.9521 

0.009 

0.9265 

0.9255 

0.013 

0.8912 

0.8899 

0.019 

0.8371 

0.8354 

0,026 

0.7553 

0.7534 

0.036 

0.6243 

0.6221 

0.049 

0.4085 

0.4049 

0.068 

0.0117 

0.0071 

0.098 

-0.5293 

-0.5345 

0.139 

-0.7352 

-0.7401 

0.185 

-0.6875 

-0.6940 

0.232 

-0.6297 

-0.6366 

0.276 

-0.5871 

-0.5930 

0.316 

-0.5537 

-0.5581 

0.353 

-0.4853 

-0.4888 

0.386 

-0.4312 

-0.4378 

0.415 

-0.4636 

-0.4682 

Longitudinal 
distance  .from 
stagnation  point 
(x/R) 

Pressure  coefficient  (Cp) 

At  side 

At  lower 
edge 

0.439 

WBEIS^M 

0.461 

0.479 

-0.4524 

-0.4448 

0.495 

-0.4462 

-0.4430 

0.515 

-0.4057 

0.569 

-0.3736 

0.676 

-0.3292 

0.835 

-0.2889 

1.041 

-0.2312 

1.289 

-0.1883 

-0.1842 

1.573 

-0.1545 

-0.1506 

1.887 

-0.1309 

-0.1263 

2.222 

-0.0827 

-o.om 

2.569 

-0.0616 

-0.0575 

2.922 

-0.0728 

-0.0694 

3.269 

-0.0717 

-0.0689 

3.604 

-0.0639 

-0.0617 

3.918 

-0.0527 

-0.0516 

4.203 

-0.0478 

-0.0481 

4.452 

-0.0487 

-0.0517 

4.659 

-0.0477 

-0.0533 

4.818 

-0.0412 

-0.0487 

4.927 

-0.0338 

-0.0458 

4.982 

-0.0231 

-0.0287 

5.002 

-0.0245 

-0.0200 

Distances  are  non-dimensional  by  maximum  head  form  radius. 
Pressure  coefficients  are  referenced  to  the  uniform  flow. 
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Table  5.  Inviscid  flow  calculation  results  for  the  20-inch  head  form  with  and  without  tunnel  walls. 


Table  5a.  Pressure  coefficients  on  the  head  form  surface. 


Longitudinal 
distance  from 
stagnation  point 
(x/R) 

Pressure  coefficient  (Cp) 

Head  form 
with  walls 

Head  form 
without 
walls 

0.000 

WKBSM 

0.000 

0.000 

0.000 

0.9965 

0.001 

0.9933 

BESB 

0.9881 

BBB 

0.9806 

0.004 

BSiB 

0.9694 

0.006 

kbb 

0.9521 

0.009 

0.9265 

0.9265 

0.013 

0.8906 

0.8910 

0.019 

0.8365 

0.8370 

0.026 

0.7550 

0.7553 

0.036 

06233 

0.6243 

0.049 

0.4063 

0.4087 

0.068 

0.0099 

0.0126 

0.098 

-0.5319 

•0.5m 

0.139 

-0.7407 

-0.7342 

0.185 

-0.6935 

0.232 

-0.6338 

0.276 

-0.5943 

-0.5874 

0.316 

-0.5586 

-0.5509 

0.353 

-0.4867 

-0.4822 

0.386 

-0.4372 

-0.4344 

0.415 

-0.4729 

-0.4648 

Longitudinal 
distance  from 
stagnation  point 
(x/R) 

Pressure  coefficient  (Cp) 

Head  form 
with  walls 

Head  form 
without 
walls 

0.439 

-0.4647 

0.461 

-0.4306 

0.479 

-0.4490 

0.495 

-0.4454 

-0.4405 

0.515 

-0.3996 

0.569 

-0.3732 

0.676 

-0.3375 

-0.3296 

0.835 

-0.2978 

-0.2891 

1.041 

-0.2398 

-0.2314 

1.289 

-0.1935 

-0.1840 

1.573 

-0.1606 

-0.1498 

1.887 

-0.1379 

-0.1263 

2.222 

-0.0883 

-0.0787 

2.569 

-0.0677 

-0.0577 

2.922 

-0.0828 

-0.0691 

3.269 

-0.0836 

-0.0684 

3.604 

-0.0772 

-0.0616 

3.918 

-0.0684 

-0.0512 

4.203 

-0.0657 

-0.0479 

4.452 

-0.0685 

-0.0515 

4.659 

-0.0701 

-0.0531 

4.818 

-0.0646 

-0.0482 

4.927 

-0.0621 

-0.0457 

4.982 

-0.0389 

-0.0253 

5.000 

-0.0363 

-0.0224 

Distances  are  non-dimensional  by  maximum  head  form  radius. 
Pressure  coefficients  are  referenced  to  the  uniform  flow. 


22 


Table  5.  (Continued). 


Table  5b.  Pressure  coefficients  along  the  strut  centerline. 


Vertical  distance 
along  strut  centerline 
(h/H) 

Pressure  coefficient  (Cp)  | 

Strut 

witb  walls 

Strut 

without  walls 

Bottom  0.00 

-0.5161 

0.10 

-0.5399 

0.20 

-0.5421 

0.30 

-0.5323 

0.40 

-0.5282 

0.50 

-0.5248 

-0.4830 

0.60 

-0.5164 

-0.4685 

0.70 

-0.5029 

-0.4454 

0.80 

-0.4944 

-0.4203 

0.90 

-0.4867 

-0.3862 

Top  1.00 

-0.4604 

-0.3248 

Vertical  distances  are  non-dimensional  by  overall  strut  height. 
Pressure  coefficients  are  referenced  to  the  uniform  flow. 
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Table  6.  Inviscid  flow  calculation  results  for  the  2-inch  and  10-inch  head  forms. 


Longitudinal 
distance  from 
stagnation  point 
(x/R) 

Pressure  coefficient  (Cp) 

2-incb  bead 
form 

10-incb  bead 
form 

0.000 

0.000 

0.000 

wEBm 

0.000 

0.001 

0.001 

mm\m 

0.003 

0.9803 

mBSM 

0.004 

0.9694 

■SI 

0.006 

0.9521 

0.009 

0.9265 

0.9265 

0.013 

0.8906 

0.8910 

0.019 

0.8365 

0.8370 

0.026 

0.7550 

0.7553 

0.036 

0.6233 

0.6243 

0.049 

0.4063 

0.4087 

0.068 

0.0099 

0.0126 

0.098 

-0.5319 

-0.5277 

0.139 

-0.7407 

-0.7342 

0.185 

-0.6935 

0.232 

-0.6338 

0.276 

-0.5943 

-0,5874 

0.316 

-0.5586 

-0.5509 

0.353 

-0.4867 

-0.4822 

0.386 

-0.4372 

-0.4344 

0.415 

-0.4729 

-0.4648 

Longitudinal 
distance  from 
stagnation  point 
(x/R) 

Pressure  coefficient  (Cp) 

2-incb  bead 
form 

10-incb  bead 
form 

0.439 

WEBBER 

-0.4647 

0.461 

-0.4306 

0.479 

-0.4542 

-0.4490 

0.495 

-0.4454 

-0.4405 

0.515 

-0.4027 

-0.3996 

0.569 

-0.3781 

-0.3732 

0.676 

-0.3375 

-0.3296 

0.835 

-0.2891 

1.041 

-0.2314 

1.289 

-0.1840 

1.573 

-0.1606 

-0.1498 

1.887 

-0.1379 

-0.1263 

2.222 

-0.0883 

-0.0787 

2.569 

-0.0677 

-0.0577 

2.922 

-0.0828 

-0.0691 

3.269 

-0.0836 

-0.0684 

3.604 

-0.0772 

-0.0616 

3.918 

-0.0684 

-0.0512 

4.203 

-0.0657 

-0.0479 

4.452 

-0.0685 

-0.0515 

4.659 

-0.0701 

-0.0531 

4.818 

-0.0646 

-0.0482 

4.927 

-0.0621 

-0.0457 

4.982 

-0.0389 

-0.0253 

5.000 

-0.0363 

-0.0224 

Distances  are  non-dimensional  by  maximum  head  form  radius. 
Pressure  coefficients  are  referenced  to  the  uniform  flow. 
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Table  7.  Viscous  flow  calculation  results  for  the  20-inch  model. 


Longitudinal 
distance  from 
stagnation  point 
(x/R) 

Pressure  coefficient  (Cp)  | 

At  side 

At  lower 
edge 

0.000 

3.2095 

4.8852 

0.001 

2.9300 

4.2880 

0.003 

2.0230 

3.8952 

0.006 

1.3683 

2.5691 

0.011 

0.8734 

1.5403 

0.017 

0.7396 

0.8620 

0.5729 

0.4874 

0.2865 

0.4268 

0.048 

-0.1212 

0.3668 

0.069 

-0.5559 

0.0101 

0.102 

-0.9352 

-0.4142 

0.151 

-1.1045 

-0.7568 

0.209 

-1.0382 

-1.0106 

0.280 

-0.9125 

-1.1480 

0.369 

-0.7634 

-1.2256 

0.477 

-0.6443 

-1.2399 

0.602 

-0.5790 

-1.0810 

0.739 

-0.5483 

-0.8974 

0.884 

-0.5260 

-0.6793 

1.028 

-0.4812 

-0.5049 

1.167 

-0.4041 

-0.3135 

1.294 

-0.2999 

-0.1855 

1.414 

-0.2072 

-0.0728 

1.539 

-0.1421 

0.0663 

1.676 

-0.1028 

0.1639 

1.823 

-0.0781 

0.2770 

Longitudinal 
distance  from 
stagnation  point 
(x/R) 

Pressure  coefficient  (Cp) 

At  side 

At  lower 
edge 

1.980 

2.146 

2.320 

0.4603 

2.500 

-0.0476 

0.4629 

2.683 

-0.0469 

0.4363 

2.869 

-0.0500 

0.4093 

3.055 

0.3616 

3.238 

0.3239 

3.416 

0.2953 

3.588 

-0.0683 

3.752 

-0.0559 

3.907 

-0.0404 

Wmm 

4.052 

-0.0383 

mSSSM 

4.185 

-0.0631 

4.308 

-0.1081 

0.2161 

4.420 

-0.1452 

0.2140 

4.521 

-0.1552 

0.2145 

4.612 

-0.1434 

0.2175 

4.693 

-0.1220 

0.2073 

4.765 

-0.1017 

0.1467 

4.828 

-0.0912 

0.0443 

4.884 

-0.0567 

-0.1060 

4.933 

0.0408 

-0.3337 

4.976 

0.1336 

-0.4517 

5.016 

0.1984 

-0.3373 

Distances  are  non-dimensional  by  maximum  head  form  radius. 
Pressure  coefficients  are  referenced  to  the  uniform  flow. 
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L _ I 

Fig.  la.  20- inch  head  form  model  and  strut 


Fig.  1.  Profile  drawings  of  the  20-inch,  10-inch,  and  2-inch  head  form  models. 
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Fig.  2.  Computational  analysis  processes  and  associated  computer  programs. 
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Fig.  3.  Shaded  picture  of  the  20-inch  model  and  strut. 


Fig.  4,  Block  structure  for  the  20-inch  model  in  the  LCC. 


Fig.  5.  Selected  block  faces  for  the  20-inch  model  in  the  LCC. 


Fig.  6.  Panel  distribution  used  for  inviscid  flow  calculations  for  the  20-inch  model  and  strut. 
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Fig.  9.  Setup  for  LDV  experiments  with  the  20-inch  model. 


Mean  longitudinal  velocity  (knots)  Mean  longitudinal 


2  20.00 


I 


18.00 


17.00 


Turbulence  =  21.5% 


32.85 


32.80- 


Number  of  samples 

Fig.  10a.  Samples  taken  near  the  tunnel  wall. 


Turbulence  =  1.67% 


32.75 


32,70- 


Mean  =  32.72  kn 


32.65- 


32.60- 


0  50  100  150  200  250  300  350  400 

Number  of  samples 

Fig.  10b.  Samples  taken  at  the  upper  edge  of  the  model,  near  the  model  surface. 


Fig.  10.  Cumulative  mean  longitudinal  velocities  versus  sample  size. 
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Number  of  samples 
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Longitudinal  velocity  (knots) 


Fig.  1  la.  Samples  taken  near  the  tunnel  wall. 


29  30  31  32  33  34  35  36 

Longitudinal  velocity  (knots) 


Fig.  lib.  Samples  taken  at  the  upper  edge  of  the  model,  near  the  model  surface. 


Fig.  11.  Distribution  of  sample  longitudinal  velocities. 
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Fig.  12.  Experimental  velocity  profiles  for  the  20-inch  model. 


Longitudinal  velocity  Longitudinal  velocity  (V^^/V^) 


0.0  0.1  0.2  0.3  0.4  0.5  0.6  0.7 

Distance  from  model  surface  (y/R) 

Fig.  12b.  Profile  around  girth,  at  lower  edge. 


0.0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8 

Distance  from  model  surface  (y/R; 

Fig.  12c.  Profile  along  upper  edge,  16  inches  aft  of  stagnation  point. 


Fig.  12.  (Continued.) 
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Pressure  coefficient  (C_)  Pressure  coefficient  (C„) 


Fig.  13a.  Pressure  coefficients  on  the  head  form  surface. 


0.0  1.0  2.0  3.0  4.0  5.0  6.0  7.0  8.0  9.0  10.0  11.0  12.0 


Longitudinal  distance  from  stagnation  point  (x/R) 

Fig.  13b.  Pressure  coefficients  on  the  length  of  the  model  surface. 


Fig.  13.  Inviscid  flow  calculation  results  for  the  20-inch  model. 
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Pressure  coefficient  (C_)  Pressure  coefficient  (C_) 
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0.0  1.0  2.0  3.0  4.0  5.0  6.0  7.0  8.0  9.0  10.0  11.0  12.0 

Longitudinal  distance  from  stagnation  point  (x/R) 

Fig.  14a.  Pressure  coefficients  on  the  model  surface. 
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—  Schiebe  results 

—  20-inch  model  without  tunnel  walls 

.  20-inch  model  with  tunnel  walls 

— 

1 

/ 

/ 

Schiebe  results  not  available  for  x/R  >  S.  12561. 

— 

Bottom  Vertical  distance  along  strut  centerline  (h/H)  Top 

of  strut  of  strut 


Fig.  1 4b.  Pressure  coefficients  on  the  strut  centerline. 


Fig.  14.  Inviscid  flow  calculation  results  for  the  20-inch  model  with  and  without  tunnel  walls. 
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Pressure  coefficient  (C_)  Pressure  coefficient 
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Fig.  17.  Viscous  and  inviscid  flow  calculation  results  for  the  20-inch  model. 
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Fig.  18.  Viscous  flow  calculation  results  for  the  20-inch  model. 
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Fig.  19.  Viscous  flow  calculation  results  versus  experimental  data. 


Distance  from  model  surface  (y/R) 

Fig.  19b.  Profile  around  girth,  at  lower  edge. 
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Fig.  19c.  Profile  along  upper  edge,  16  inches  aft  of  stagnation  point. 


Fig.  19.  (Continued). 


41 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK. 


42 


APPENDIX  A 

DETAILED  PROCEDURE  FOR  COMPUTATIONAL  ANALYSIS 


43 


GEOMETRIC  MODELLING 

The  first  step  of  the  computational  analysis  is  to  develop  numerical  models.  Offsets  are 
gathered  from  any  available  source;  lines  plans,  model  shop  drawings,  or  electronic  offsets. 
Typically,  commercial  software  packages  for  personal  computers,  such  as  MACS  UR  F  or 
FASTSHIP,  are  used  to  generate  the  three-dimensional  numerical  surfaces  representing  the 
desired  physical  model,  using  conventional  Bsplines  or  non-uniform,  rational  B-splines 
(NURBs).  Conventional  B-splines  cannot  precisely  represent  the  shape  of  conic  sections,  but 
NURBs  have  this  capability. 

Sometimes,  in  a  case  such  as  the  cavitation  head  forms,  the  models  consist  primarily  of 
known  geometrical  shapes  (e.g.,  nearly  flat  fins  on  a  body  of  revolution  and  a  tapered  ogive 
strut).  It  is  convenient  to  generate  the  models  directly  by  writing  a  simple  computer  program  to 
handle  each  specific  geometry. 

As  the  surfaces  are  defined,  a  grid  is  distributed  over  the  surface.  This  step  is  handled  by 
a  program  developed  in-house  called  REMESH.  In  instances  where  the  surfaces  are  defined 
directly,  the  surface  is  represented  as  a  grid  from  the  beginning,  and  this  process  is  not  an 
individual  step.  The  intersections  of  the  various  surfaces  are  computed  and  the  surfaces  are 
trimmed  to  meet  at  the  intersections  without  overlapping.  Using  another  computer  program 
(NWIRE,  also  developed  in-house),  shaded  pictures  are  generated  and  visually  inspected  to 
ensure  accurate  surface  generation. 

Different  types  of  numerical  grids  can  be  generated  around  the  numerical  models,  for 
different  methods  of  computational  analysis,  generally  viscous  flow  or  inviscid  flow  calculations. 
For  the  viscous  flow  analysis,  a  volume  grid  is  developed,  typically  extending  from  the  model 
surface  to  some  outer  boundary.  For  the  inviscid  flow  analysis,  a  panelled  surface  is  generated 
by  redistributing  the  points  on  the  surface  grid. 
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VISCOUS  FLOW  CALCULATIONS 
Grid  Generation 

Volume  grids  for  viscous  flow  calculations  can  be  developed  using  a  computer  package 
specifically  developed  for  three-dimensional,  multi-block  grid  generation.  This  package,  called 
GRIDGEN,  was  developed  by  General  Dynamics  for  the  U.S.  Air  Force.  Using  GRIDGEN,  the 
topology  of  the  blocks  and  the  connectivity  of  the  blocks  is  defined.  Tools  in  the  GRIDGEN 
package  allow  for  control  of  the  grid  size  and  distribution  on  each  face  of  each  block. 

Once  the  face  grids  are  defined,  a  program  developed  in-house  called  TRANS3D  is  used. 
TRANS3D  determines  the  coordinates  of  the  points  on  the  interior  of  each  block  using  a  three- 
dimensional,  algebraic,  transfinite  interpolation  scheme.  Essentially,  the  grids  on  the  opposite 
faces  of  each  block  are  projected  through  the  interior  of  the  block.  By  projecting  through  in  three 
directions  (three  pairs  of  opposing  faces  on  a  six-sided  block),  points  are  located  in  space. 
These  points  outline  cells,  or  three-dimensional  volumes.  The  program  checks  each  cell  to 
ensure  that  the  cells  are  oriented  consistently  with  a  right-handed  coordinate  system.  Iterations 
of  grid  generation  are  made  to  achieve  an  acceptable  three-dimensional  volume  grid  with  a 
relatively  smooth,  and  properly  oriented,  grid  distribution. 

Solution  Technique 

The  viscous  flow  calculations  are  made  using  a  program  called  JSFLOW.  The  three- 
dimensional  volume  grid  is  developed  for  input  to  the  program,  and  ISFLOW  calculates  a  steady- 
state,  incompressible  flow  solution  using  the  three-dimensional  Reynolds-averaged  Navier- 
Stokes  (RANS)  equations.  The  three-dimensional  velocity  components  are  computed  at  each 
point  in  the  grid  through  an  iterative  process  for  the  steady-state  solution.  A  more  complete 
discussion  of  the  program  operation  is  given  in  Lin  et  al.^ 

A  grid  for  an  unappended  surface  ship  might  require  about  400,000  points,  versus  over  a 
million  for  a  fully  appended  model.  Grid  points  are  typically  concentrated  in  areas  of  high 
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Curvature  on  the  model,  or  in  areas  of  rapidly  changing  flow  conditions,  such  as  near  a  stagnation 
point.  The  head  form  models  therefore  required  about  600,000  grid  points,  due  to  the  extremely 
rapid  variation  of  flow  characteristics  near  the  forward  end  of  the  models  (e.g.,  near  the  high 
curvature  of  the  head  form  and  near  the  stagnation  point).  A  grid  of  this  size  requires  about  two- 
and-a-half  minutes  of  computational  time  per  iteration  of  the  RANS  equations  on  a  powerful 
workstation,  or  about  one  minute  on  a  super-computer.  Typically,  2500  to  3000  iterations,  or 
time-steps,  are  required  before  the  program  reaches  a  steady  solution.  A  solution,  therefore, 
requires  anywhere  from  50  to  125  hours  of  computer  time  to  be  calculated.  Since  the  grid  must  be 
refined  based  on  computational  results,  several  incomplete  runs  of  500  or  1000  time-steps  are 
generally  made  while  the  final  grid  is  developed  and  refined.  A  single  grid  can  require  a  total  of 
250  to  350  hours  of  computer  time  to  reach  a  final  solution. 

INVISCID  FLOW  CALCULATIONS 
Model  Generation 

Inviscid  flow  calculations  typically  require  a  panelled  surface  as  input,  which  is  generated 
by  redistributing  the  points  on  the  previously  generated  surface  grid.  Programs  developed  in- 
house  called  REMESH  and  REPAN  allow  control  of  the  panel  distribution,  while  ensuring  that 
the  panels  conform  to  the  surface  shape.  The  panel  distribution  must  be  adjusted  to  concentrate 
more  panels  in  areas  of  high  curvature,  radical  change,  or  areas  of  special  interest.  At  the  same 
time,  it  is  desirable  to  minimize  the  number  of  panels  in  order  to  minimize  computation  time. 
Rapid  changes  in  panel  size  or  aspect  ratio  can  lead  to  poor  computational  results. 

Solution  Technique 

The  inviscid  flow  calculations  can  be  made  using  several  different  computer  programs.  In 
this  project,  the  calculations  were  made  using  a  program  called  VFLOW.  VFLOW  does  not 
account  for  a  free  surface,  but  a  free-surface  representation  was  not  needed  to  examine  the  flow 
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in  a  closed  tunnel  such  as  the  LCC. 

A  panelled  surface  is  supplied  as  input  to  the  program,  and  the  program  then  distributes 
doublet  and  source  singularities  over  the  surface.  VFLOW  assumes  incompressible  flow.  Local 
velocity  in  three  components  and  pressure  coefficient  are  calculated  for  each  panel.  As  indicated 
by  the  name,  inviscid  flow  calculations  do  not  account  for  viscous  effects,  such  as  turbulence  and 
separation. 

An  inviscid  model  typically  has  1500  to  2500  panels,  versus  hundreds  of  thousands  to 
over  a  million  grid  points  for  a  viscous  grid.  The  solution  time  is  correspondingly  quicker.  An 
inviscid  flow  solution  can  take  anywhere  from  2  to  5  minutes  of  computational  time.  Even 
multiple  iterations  while  refining  the  model  do  not  involve  much  computer  time. 

PRESENTATION  OF  RESULTS 

The  output  produced  by  IS  FLOW  includes  a  flow  solution  file  that  contains  pressure  and 
velocities  in  three  components  for  each  point.  This  solution  file  and  the  grid  coordinate  file  were 
used  to  display  the  results  in  various  graphical  formats.  A  program  called  PLOT3D,  developed 
by  NASA  Ames  Research  Center,  is  often  used  to  produce  three-dimensional  plots  of  geometry, 
grids,  and  computational  results  as  streamlines,  vectors,  contours,  etc. 

PLOTS D  was  developed  for  compressible  flow.  Pressure  and  velocity  are  converted  to 
density,  stagnation  energy,  and  momentum  values.  PLOT3D  can  then  be  used  to  represent  these 
values  as  streamlines,  vectors,  contours,  etc.  of  various  functions.  The  PLOT3D  user's  manual^ 
contains  a  complete  list  of  available  functions. 

No  programs  were  available  specifically  for  viewing  the  inviscid  flow  calculation  results. 
However,  PLOT3D  is  versatile  enough  to  accept  a  modified  input  file  for  basic  plotting.  Plotting 
was  limited  to  the  specifically  generated  output  from  VFLOW.  No  plots  could  be  made  where  no 
information  was  present  (e.g.,  anywhere  off  the  model  or  strut  surface.) 

In  add'tion,  the  results  can  be  graphed  in  two-dimensional  formats  using  various 
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Commercially  available  graphing/plotting  packages. 
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APPENDIX  B 

EXPERIMENTAL  RESULTS  FOR  THE  20-INCH  HEAD  FORM  IN  THE  LCC 
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Table  B-1.  Experimental  data. 
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Table  B-ld.  Profile  around  girth,  at  upper  edge,  22.9  knots. 
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This  point  is  the  low  turbulence  data  point  shown  in  Figures  9b  and  10b. 
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These  profiles  are  graphed  in  Figures  11b  and  18b. 
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This  profile  is  graphed  in  Figures  1  Ic  and  18c. 


Table  B-lp.  Profile  at  centerline,  from  model  to  tunnel  wall,  14.0  knots. 
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This  profile  and  the  similar  profiles  at  different  speeds  are  graphed  in  Figure  1  la. 
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*  This  profile  and  the  similar  profiles  at  different  speeds  are  graphed  in  Figure  11a.  This  profile  is  also  graphed  in 
Figure  18a. 

t  This  point  is  the  high  turbulence  data  point  shown  in  Figures  9a  and  10a. 
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This  profile  and  the  similar  profiles  at  different  speeds  are  graphed  in  Figure  1  la. 
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